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La presente invention concerne une methode et un systeme pour estimer en 
temps reel le mode d'ecoulement d'une veine fluide polyphasique en tous points 
d'une conduite comprenant l'utilisation des reseaux neuronaux. 

Etat de la technique 

L'acheminement des hydrocarbures depuis les sites de production jusqu'aux 
unites de traitement constitue un maillon important de la chaine petroliere. G'est un 
maillon delicat en raison de la complexite des interactions entre les phases 
constituant les effluents transportes. Les operateurs ont pour objectif premier 
d'atteindre une productivity optimale dans les meilleures conditions de securite. lis 
doivent done gerer au mieux la vitesse et la temperature, pour eviter des pertes de 
charges superflues, des depots indesirables et des irregularites d'ecoulement. La 
methode generalement utilisee consiste a modeliser au mieux le transport de flux 
polyphasiques complexes de fagon a fournir a chaque instant une image des 
ecoulements dans les differentes parties de la chaine de production, tenant compte de 
la constitution precise de l'effluent, les debits et pressions et les modes d'ecoulement. 

D existe actuellement differents outils logiciels de simulation du transport de 
flux polyphasiques complexes, permettant a une stade precoce de concevoir des 
equipements de production adaptes. 

Par les brevets US 5 550 761, FR 2.756.044 (US 6_028 992} et FR 2 756 045 
(US 5 960 187) du demandeur, notamment, on connait des methodes et outils de 
modelisation permettant de simuler le transport de flux polyphasiques complexes en 
regime permanent ou transitoire et capables de prendre en compte des phenomenes 
d'instabilite qui se produisent du fait de la geometrie irrdguliere du terrain ou passe la 
conduite ou de la topographie de celle-ci, que les specialistes designent par « terrain 
slugging » ou « severe slugging ». 

La complexite des outils de simulation est a 1' image de celle des phenomenes 
modelises. Precision et performances ne peuvent etre obtenues qu'apres un temps de 
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modelisation relativement important qui s'avere difficilement compatible avec une 
gestion en temps reel. C'est pourquoi ces outils de modelisation ne peuvent en Tetat 
servir a la gestion en temps reel de la production. Dans cette optique, il est ainsi 
apparu necessaire de mettre en ceuvre des methodes de modelisation permettant un 
5 bon compromis entre rapidite de calcul et precision des resultats. 

L' invention a pour objet essentiel une methode permettant seule ou en 
parallele avec les methodes de modelisation ci-dessus, de gerer en temps reel des 
parametres d'une circulation de fluides en utilisant des reseaux neuronaux. 

Les reseaux neuronaux definissent on le rappelle un mode de traitement de 

10 donnees simulant le fonctionnement des systemes de neurones biologiques. Dans de 
tels reseaux, un element realise un calcul relativement simple tel qu'une somme 
ponderee des signaux presents a ses entrees appliquees a une fonction non-lineaire, 
qui determine l'etat de sa sortie. On utilise un nombre important de tels elements 
interconnects en serie et en parallele. Un choix convenable des facteurs de 

15 ponderation permet au reseau de realiser des fonctions complexes. Les reseaux dits a 
retropropagation par exemple utilisent des couches multiples d'elements definis 
ci-dessus. L' adaptation d'un tel reseau a une tache precise est faite en «entrainant » 
le reseau sur un certain nombre d'exemples et en ajustant les facteurs de ponderation 
pour chaque element aux valeurs qui conviennent. On presente des valeurs d' entree 

20 au reseau, on analyse la valeur de sortie produite par le reseau et on modifie les 
facteurs de ponderation pour minimiser au mieux l'ecart entre la valeur effective a la 
sortie et la valeur attendue dans l'exemple choisi. Apres un entrainement suffisant, le 
reseau est adapte a repondre a de nouvelles valeurs d'entree pour lesquelles la valeur 
de sortie n'est pas connue a priori et a produire une valeur de sortie adaptee. Dans 

25 son principe, un reseau de neurones procede selon une methode de regression non 
lineaire, d'autant plus performante par rapport aux methodes classiques. 

On utilise de tels reseaux dans de nombreux domaines tels que la 
reconnaissance d'images, la resolution de problemes d'optimisation, etc. 
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Initialement, le reseau de neurones est une methode adaptee k la classification 
automatique, d'ou son utilisation en particulier pour la reconnaissance de formes. 
Pour ces applications, deux types de reseaux peuvent etre mis en oeuvre, les reseaux 
MLP (Multi Layer Perceptron) ou les reseaux de Kohonen, bien connus des 
5 specialistes. 

L'art anterieur dans le domaine des reseaux neuronaux est illustre par les 
references suivantes : 

Dreyfus G., «Les reseaux de neurones » ; Mecanique Industrielle , et materiaux, 
n°51, sept. 98 ; 

10 - Lippman R. P., An Introduction to Computing with Neural Nets ; IEEE ASSP 
Magazine, April 1987 ; or 

Pinkus A., Approximation Theory of the MLP Model in Neural Networks ; Acta 
Numerica 1999. 

Les reseaux sont maintenant egalement employes pour la modelisation non- 
15 lineaire de donnees statiques ou de processus dynamiques. Dans ce cas, ce sont 
majoritairement les reseaux MLP qui sont utilises. Cette approche concerne 
actuellement des domaines d' applications comme, par exemple, la detection 
d' anomalies ou la prediction boursiere. 

Un exemple d' utilisation de reseaux neuronaux est decrit par exemple dans le 
20 brevet FR -A- 2.786.568 du demandeur. 

La methode selon Finvention 

Dans le contexte de la gestion de production, les phenomenes a modeliser en 
temps reel sont fortement non lineaires et les parametres intervenant dans la 
modelisation fort nombreux. 



4 



L' invention a pour objet une methode pour estimer en temps reel le mode 
d'ecoulement en tous points d'une conduite de structure que Ton peut definir par un 
certain nombre de parametres de structure, d'une veine de fluide polyphasique 
definie par plusieurs grandeurs physiques et comprenant au moins une phase liquide 
5 et au moins une phase gazeuse, caracterisee en ce que la modelisation du mode 
d'ecoulement comporte : 

la formation d'un reseau neuronal non lineaire avec une couche d'entree avec 
autant d'entrees que de parametres de structure et de grandeurs physiques, une 
couche de sortie avec autant de sorties que de grandeurs necessaires a 
10 Testimation du mode d'ecoulement et au moins une couche intermediaire ; 

la constitution d'une base d'apprentissage avec des tables predefinies reliant 
differentes valeurs obtenues pour les donnees de sortie aux valeurs 
correspondantes des donnees d'entree ; et 

la determination par iterations des facteurs de ponderation de la fonction 
15 d' activation permettant de relier correctement les valeurs dans les tables des 

donnees d'entree et de sortie. 

La methode comporte en outre de preference une analyse des donnees de 
sortie du reseau de neurones permettant de trier, parmi les valeurs des donnees de 
sortie du reseau de neurones, les seules donnees pertinentes a prendre en compte dans 
20 la determination iterative des coefficients de ponderation de la fonction d' activation. 

Un reseau neuronal traduit une fonction mathematique continue entfe des 
donnees d'entree et de sortie. Comme il produit en continu des resultats par des 
calculs directs reliant en permanence des donnees d'entree et de sortie, on evite par 
consequent les discontinuity liees a la resolution de systemes d'equations, d'ou un 
25 gain de temps appreciable et une simplification des resultats. Le caractere non 
lineaire des fonctions d' activation des reseaux de neurones les rend parfaitement 
adaptables a une modelisation de systemes non lineaires. 



La methode comporte par exemple la formation d'un reseau totalement 
connecte et 1' utilisation de neurones de sortie lineaires. 

Comme fonction d' activation, on choisit par exemple une fonction identite. 

Le systeme selon 1' invention permet d'estimer en temps reel le mode 
d'ecoulement en tous points d'une conduite de structure que Ton peut definir par un 
certain nombre de parametres de structure, d'une veine de fluide polyphasique 
definie par plusieurs grandeurs physiques et comprenant au moins une phase liquide 
et au moins une phase gazeuse. II comporte : 

des moyens pour determiner les caracteristiques d'un reseau neuronal non lineaire 
avec une couche d'entree avec autant d'entrees que de parametres de structure et 
de grandeurs physiques, une couche de sortie avec autant de sorties que de 
grandeurs necessaires a l'estimation du mode d'ecoulement et au moins une 
couche intermediaire ; 

des moyens de stockage d'une base d'apprentissage avec des tables predefinies 
reliant differentes valeurs obtenues pour les donnees de sortie aux valeurs 
correspondantes des donnees d' entree ; et 

des moyens pour determiner par iterations des facteurs de ponderation d'une 
fonction d' activation permettant de relier correctement les valeurs dans les tables 
des donnees d' entree et des donnees de sortie. 

Le systeme comporte de preference des moyens d' analyse des donnees de 
sortie du reseau de neurones permettant de trier, parmi les valeurs des donnees de 
sortie du reseau de neurones, les seules donnees pertinentes a prendre en compte dans 
la determination iterative des coefficients de ponderation de la fonction d' activation. 

Le systeme developpe par cette technique peut etre utilise par exemple en 
remplacement d'un outil developpe par une technique classique (et done relativement 
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lent) dans le but de fournir des resultats compatibles avec une gestion des flux en 
temps reel. Les termes de cette mise en ceuvre sont developpes ci-dessous. 

Presentation des figures 

La methode sera mieux comprise a la lecture de la description ci-apres en se 
5 referant aux dessins sur lesquels : 

la Fig.l qui montre un exemple de reseau neuronal construit pour relier des 
donnees hydrodynamiques d' entree a des donnees de sortie portant sur des modes 
d'ecoulement ; et 

les Fig. 2 et 3 illustrent deux cas de voisinages topologiques differents conduisant 
10 a une qualification differente de donnees de sortie dans un processus de 

validation du reseau neuronal. 

Description detaillee 

1) Contexte 

On considere une circulation de fluides polyphasiques dans une conduite avec 
15 au moins une phase liquide et au moins une phase gazeuse et Ton cherche a 
construire un reseau neuronal permettant, a partir d'un certain nombre de donnees 
d' entree geometriques et physiques relatives a la conduite et de donnees physiques 
sur les fluides, de donner a chaque instant et pour chaque section de la veine fluide 
une estimation du mode d'ecoulement. 

20 2) entrees et sorties 

Les donnees d'entree sont par exemple : 

des donnees geometriques du reseau : diametre, rugosite et angle d'inclinaison du 
conduit ; 
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des donnees d'entree qualifiant les effluents : masse volumique de gaz, masse 
volumique de liquide, viscosite du gaz, viscosite du liquide, etc. ; 

des donnees d'entrees caracterisant le melange : tension superficielle gaz/liquide, 
fraction volumique de gaz, vitesse barycentrique du melange. 

5 Les donnees que le reseau va evaluer et delivrer sur deux sorties principales 

sont : 

dV, la difference de vitesse entre gaz et liquide ; et 

P, la fraction d'ecoulement stratifie dans la partie de conduit dont on souhaite 
determiner le type d'ecoulement ; (3 e [0 ; l] . 

10 D'autres grandeurs qualifiant le type d'ecoulement peuvent etre calculees a 

partir de ces deux sorties. 

3) structure du reseau 

Pour relier entre elles les donnees d'entree et de sortie, on construit un reseau 
de neurones de preference de type MLP bien connu des gens de Fart, car il est 

15 particulierement bien adapte a la modelisation de phenomenes physiques. En effet, sa 
structure permet de decrire tout a la fois les composantes dynamiques et statiques des 
phenomenes, en fixant meme si necessaire certains de ses parametres a une valeur 
reifiee, et done physiquement representative. Ainsi, dans l'exemple, la connaissance 
d'equations physiques regissant les ecoulements permet d'enrichir le reseau et de 

20 l'adapter au mieux aux phenomenes physiques qu'il modelise. 

Le reseau de neurones comporte (Fig.l) trois couches : la couche d'entree 
composee de dix neurones correspondant aux dix donnees (citees ci-dessus) du 
modele physique complet, une couche de sortie de deux neurones correspondant aux 
deux parametres recherches dV et (3, et une couche intermediaire, dite couche cachee, 
25 dont le nombre de neurones N c est optimise. Le reseau est totalement connecte. La 
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non linearite de ce reseau est obtenue par une fonction d' activation sigmoide 
regissant le comportement des neurones de la couche cachee. Les neurones de la 
couche de sortie peuvent etre choisis lineaires ou non-lineaires. La fonction 
d' activation peut etre par exemple la fonction identite. 

5 4) Apprentissage : principe et mise en oeuvre dans l'exemple 

a) Principe 

Les poids de cette structure sont determines a Tissue d'une phase 
d'apprentissage ; cette phase consiste a nourrir le reseau d'un ensemble de donnees 
constituant la base d' apprentissage du reseau, et a optimiser les poids du reseau par 
10 minimisation des erreurs constatees pour l'ensemble des echantillons de la base, 
entre les donnees de sortie issues du calcul du reseau et les donnees attendues a la 
sortie, donnees par la base. Les erreurs peuvent etre les erreurs absolues entre les 
grandeurs d'entree et de sortie ou les erreurs relatives, selon la performance desiree 
pour le reseau. 

15 Les facultes de generalisation du reseau sont ensuite testees sur sa capacite a 

bien calculer les deux sorties pour des entrees qui lui sont inconnues. 

b) Mise en oeuvre 

En pratique, outre la difficulte du choix des differents elements constituant le 
reseau, la mise en place du reseau passe par une analyse poussee des donnees 

20 constituant la base d'apprentissage. En effet, meme si un reseau est convenablement 
elabore pour un probleme donne, il peut donner in fine de mauvais resultats dans la 
mesure ou la base d'apprentissage ; ; qui lui est fournie contient des elements qui 
perturbent I'optiniisation de l'ensemble de ses poids. Dans le cas particulier des 
ecoulements de fluides, le probleme est crucial ; de fait, les calculs de modes 

25 d'ecoulements realises a partir de modeles complets et precis generent une base de 
donnees fortement inhomogene : par exemple, un resultat stipulant un ecoulement 
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stratifie peut etre « noye »' dans un voisinage de points representant un ecoulement 
disperse : ce point correspond done soit un phenomene physique que le r£seau de 
neurones doit pouvoir retrouver, soit un « singleton » que le reseau doit au contraire 
lisser pour assurer la continuite du rtsultat. La difficulte de F analyse prealable des 
5 donndes reside done dans la determination de la nature de tels points, et, a plus 
grande echelle, des differents points de la base. 

Pour cela, une methode d' analyse a posteriori des resultats du reseau a ete 
elaboree, qui identifie et distingue les donnees disponibles representant les 
phenomenes physiques particuliers a modeliser, des resultats de calculs a ignorer. 

10 Cette methode consiste a fabriquer des « compacts » autour de chacun des points 
singuliers consideres dans l'espace a n dimensions des entrees (n = 10 dans 
l'exemple decrit), afin d'evaluer le degre d'isolement de leur comportement ; ainsi, 
pour chaque « boule » creee sont comptabilises les points singuliers et les points non 
singuliers, la proportion des premiers dans le compact donnant alors une mesure de 

15 leur isolement. Les termes « compacts et « boule » sont ici pris au sens topologique 
du terme. 

Ainsi, le point SI de la Fig.2 qui est au centre d'une boule et separe du point 
singulier S2 le plus proche par un voisinage de nombreux points reguliers Mi, peut 
etre considere comme un singleton et done etre regularise. Le point SI de la Fig.3, au 
20 contraire, a pour voisinage immediat un ou plusieurs points egalement singuliers S2, 
S3,... Sk. Dans un tel cas, il faut considerer que les points SI a Sk sont la traduction 
d'un phenomene physique particulier dont il faut tenir compte par une adaptation du 
reseau neuronal. 

Cette methode de selection permet d' identifier les points necessaires pour 
25 informer suffisamment le reseau lors de la phase d'apprentissage. La base de points 
adequats ainsi constitute est alors capable de le renseigner afin d'optimiser au mieux 
ses parametres. 
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5) Resultats 

Avec les outils et methodes decrits ci-dessus, on obtient : 

une base de donnees nettoyee de tout point a comportement anormal et contenant 
uniquement les points riches en information pour cerner au mieux le systeme 
physique a identifier ; 

un reseau capable d'interpoler les differents types d'ecoulements rencontres, afin 
de ne plus generer de discontinuity, et egalement apte a tenir compte des 
phenomenes physiques particuliers dont on souhaite la modelisation malgre leur 
comportement atypique ; et finalement 

un outil donnant en temps reel une estimation des principales informations 
d' hydrodynamique. 
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REVENDICATIONS 



1) Methode pour estimer en temps reel le mode d'ecoulement, en tous points 
d'une conduite de structure que Ton peut definir par un certain nombre de parametres 
de structure, d'une veine de fluide polyphasique definie par plusieurs grandeurs 
physiques et comprenant au moins une phase liquide et au moins une phase gazeuse, 
caracterisee en ce qu'elle comporte : 

la formation d'un reseau neuronal non lineaire avec une couche d'entree avec 
autant d'entrees que de parametres de structure et de grandeurs physiques, une 
couche de sortie avec autant de sorties que de grandeurs necessaires a 
l'estimation du mode d'ecoulement et au moins une couche intermediaire ; 

la constitution d'une base d'apprentissage avec des tables predefinies reliant 
differentes valeurs obtenues pour les donnees de sortie aux valeurs 
correspondantes des donnees d'entree ; et 

la determination par iterations des facteurs de ponderation de la fonction 
d' activation permettant de relier correctement les valeurs dans les tables des 
donnees d'entree et de sortie. 

2) Methode selon la revendication 1, caracterisee en ce qu'elle comporte en 
outre une analyse des donnees de sortie du reseau de neurones permettant de trier, 
parmi les valeurs des donnees de sortie du reseau de neurones, les seules donnees 
pertinentes a prendre en compte dans la determination iterative des coefficients de 
ponderation de la fonction d' activation. 

3) Methode selon la revendication 1 ou 2, caracterisee en ce que Ton forme 
un reseau totalement connecte. 
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4) Methode selon Tune des revendications 1 a 3, caracterisee en ce que les 
neurones de sortie sont lineaires. 

5) Methode selon Tune des revendications 1 a 3, caracterisee en ce que la 
fonction d'activation est une fonction identite. 

5 6) Systeme pour r estimation en temps reel du mode d'ecoulement, en tous 

points d'une conduite de structure que Ton peut definir par un certain nombre de 
parametres de structure, d'une veine de fluide polyphasique d6finie par plusieurs 
grandeurs physiques et comprenant au moins une phase liquide et au moins une 
phase gazeuse, caracterise en ce qu'il comporte : 

10 - des moyens pour determiner les caracteristiques d'un reseau neuronal non lineaire 
avec une couche d'entree avec autant d' entrees que de parametres de structure et 
de grandeurs physiques, une couche de sortie avec autant de sorties que de 
grandeurs necessaires a Testimation du mode d'ecoulement et au moins une 
couche intermediate ; 

15 - des moyens de stockage d'une base d'apprentissage avec des tables predefinies 
reliant differentes valeurs obtenues pour les donnees de sortie aux valeurs 
correspondantes des donnees d'entree ; et 

des moyens pour determiner par iterations des facteurs de ponderation d'une 
fonction d' activation permettant de relier correctement les valeurs dans les tables 
20 des donnees d'entree et des donnees de sortie. 

7) Systeme selon la revendication 6, caracterise en ce qu'il comporte des 
moyens d' analyse des donnees de sortie du reseau de neurones permettant de trier, 
parmi les valeurs des donnees de sortie du reseau de neurones, les seules donnees 
pertinentes a prendre en compte dans la determination iterative des coefficients de 
25 ponderation de la fonction d' activation. 
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